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Abstract 

We study a model for pair interaction U of DNA molecules generated by the discrete dipole 
moments of base-pairs and the charges of phosphate groups, and find noncommutative group of 
eighth order S of symmetries that leave U invariant. We classify the minima using group S and 
employ numerical methods for finding them. The minima may correspond to several cholesteric 
phases, as well as phases formed by cross-like conformations of molecules at an angle close to 90° , 
"snowflake phase". The results depend on the effective charge Q of the phosphate group which can 
be modified by the polycations or the ions of metals. The snowflake phase could exist for Q above 
the threshold Qc- Below Qc there could be several cholesteric phases. Close to Qc the snowflake 
phase could change into the cholesteric one at constant distance between adjacent molecules. 
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I. INTRODUCTION. 



The DNA cholesterics have drawn considerable attention as unusual states of condensed 
matter, . These liquid-crystalline phases are of interest from various points of view, 

[2I. For one thing there is a certain affinity, as packaging is concerned, between them and 
viruses and chromosomes, for another they could be useful for constructing nano-devices. 



for example sensors. 



Cholesteric assemblies of the DNA have been observed in vivo and extensively studied in 



vitro, 



1^. The experimental methods for obtaining molecules of the DNA in condensed 
state employ either solutions of polycations, or of polymers separating phases containing 
molecules of the DNA. Livolant et al, [ol, used polycations and found the spatial 

characteristics of the location of DNA molecules, [ll|. Yevdokimov et al, [2], have used 
the PEG pellets for the formation of liquid crystalline phases of the DNA, and found the 
polymorphism by observing a change in the sign of cholesteric spiral. 

By now it is generally accepted that the liquid crystalline phases of the DNA are essen- 



tional properties of the DNA playing the crucial role, [13|] 



tially of the electrostatic nature due to the helical distribution of charge, [1^, the conforma 



The concept of the charged 



duplex has made for understanding the physics of the phases, and enabled to explain t 



ne 



emergence of the cholesteric structure and its relation to the ion strength of solutions, 151] • 



But the existence of different cholesteric phases of the DNA is still a challenge to the 
current theory. The theoretical study of pair-interaction U caused by charges of phosphate 
groups, dipoles of base-pairs, and ions of ambient solution, can make for understanding the 
phenomenon. Electric charges of molecules of the DNA are of two kinds: those due to the 
phosphate groups and those of the base-pairs, which are generally described by dipole mo- 
ments, jl^ • Kim, , used the helix symmetry of the set of dipole moments, for describing 
a torque that may deform the parallel conformation of two adjacent molecules and result 
in a cholesteric phase of the system (see also [IS]). The idea turned out to be quite useful 
for the physics of the DNA. The main point about its implementation is that the discrete 
structure of charges presents serious difficulties for analytical treatment. 
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II. PAIR INTERACTION OF DNA MOLECULES. 

n 

In this paper we study the multistabihty of cholesteric phases, |8i], within the framework 
of the approach indicated above, taking into account only the electrostatic forces between 
molecules of the DNA in solution. We use the discrete approximation for the electrostatic 
interaction between pairs of molecules of the DNA, [17], use extensively computer simulation, 
and employ symmetry arguments. Our model is essentially that of one dimensional lattice of 
charges and dipoles with an elementary cell of size 3.4 A that mimics the spatial conformation 
of the charges of phosphate groups and the dipoles of base-pairs. We consider short segments 
of the DNA, approximately 500 A, that is of the size of persistence length, so that to a good 
approximation one may consider them as segments of straight lines. We shall assume that 
both molecules have the same number of base-pairs that can be visualized as points on a 
straight line parallel to the axis of the molecule, one base-pair being located at the center of 
a corresponding molecule. The centers of the straight lines describing the molecules belong 
to the straight line perpendicular to the plane x-y parallel to either of them. We shall denote 
by ^ the angle between the straight lines describing the molecules. We assume that both 
molecules are of the same helicity. The latter is determined by the rotation of the frame of 
the dipole moments. Thus, we model a molecule of the DNA on a one-dimensional lattice 
having at its sites either vectors of dipoles of the base pairs or scalars of the phosphate 
charges. The case of total neutralization of phosphate charges was considered in paper jl9 |. 

The dipoles are suggested to have the helix symmetry with vr/S rotation / bp, corre- 
sponding to the structure of the ideal double helix of the DNA. The distance between the 
centers n, which is fixed, is an important parameter of the model. In what follows we use 
the distance between adjacent base-pairs, that is 3.4 A, as the unit of length, and take a 
unit of charge for which the dipole moment of 1 Debye equals 1. We perform calculations 
in the dimensionless units generated by these quantities. 

The energy of electrostatic interaction of two molecules can be cast in the sum 

ell = Uo + Udd + Udc + Ucd + Ucc (1) 

in which e is the dielectric permeability of solvent and Uq is the self energy of the pair, which 
does not influence its conformation, the first term describes the interaction between dipoles 
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of the first molecule and those of the second; the second term - dipoles of the first and 
phosphate charges of the second; the third - charges of the first and dipoles of the second; 
the fourth - charges of the first and the second. The interactions are given by the equations 



Udd(r, r 



Udc{r,r 



-up 



9{P) — {P ■ P ) - 3/i(p) 



-up 



k{p)Q 



Ucd{r',r) = e ''''k{p)Q 
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in which u is the inverse Debye length u = A~^, and 

p = |f^— f'\ 

We shall take the screening functions k{p),g{p), h{p) in Schwinger's form 

k = g = l + i'p, h = 1 -\- u p -\ — z/^p^ 
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(2) 
(3) 
(4) 
(5) 



The important point about the electrostatic interaction between molecules of the DNA 
is a wise choice of the screening factor. The common practice is to employ the Debye- 
Hiickel theory, or its modifications that might accommodate the dipole charges, the so-called 
Schwinger screening. The full treatment of this problem requires a separate investigation. 
In this paper we confine ourselves to the Debye-Hiickel and the Schwinger theories, |20l ]. 



III. SYMMETRY GROUP FOR PAIR-INTERACTION. 

In what follows we shall use the standard technic of group theory, which enables us to 
formulate the properties of minima for potential ?7 in a concise and clear form. 

It is worth noting that the pair potential U is invariant as regards the change of helicity: 
if we change the sign of angle ^ between the axes of the two molecules, at the same time as 
the sign of helicity, the potential U remains the same. Thus, there is the invariance rule as 
regards the transformation 

(right hehcity, ^) <^ (left hehcity, - ^) (6) 
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There are symmetry rules for the hehxes of the same kind. One may convince oneself, for 
example by writing explicit expressions for the interactions given by equations ([S])-© with 
the help of angles 0i,02 between the dipole moments at the centers of the two molecules 
and axis 2;, and the angle ^ introduced above, that the following transformations 

h : (01, 02, ^ (01, vr-02, e + vr) (7) 
t2 : (01, 02, ^ (vr-01, 02, e + vr) (8) 
h : (01, 02, (02 + vr, 0i + 7r, (9) 

leave the potential U invariant. The angles are defined within limits 

— Vr < 01 < TT, — TT < 02 < VT, —71 < ^ < 

values ±71 corresponding to the same configurations of the molecules. 
The transformations given by equations dTJlQ]) verify the equations 

2 2 2 "J 
tl — ^2 ~ ^3 ~ ^2 ^3 = ^3 ^1, ti ^2 = ^2 ^1, 

where id is a transformation that leaves all 0i, 02, ^ invariant. Using the above equations one 
can easily convince oneself that ^1,^2, ^3 generate a non- commutative group of 8-th order, S. 
Its maximal subgroup ?i is a normal subgroup of 4-th order, commutative, and generated 
by the transformations 

fl=h, /2 = tlt2t3 (10) 

Elements /i, /2 in its turn generate subgroups Hi and 7^2 of TC, respectively. It is worth 
noting that Hi, H2 are of second order, both. They are conjugate subgroups of S, that is 
for an element g of S we have fi = g^^ /2 g, or we may state Hi = g^^ H2 g, in the notations 
of group theory, which can be cast in the form of the diagram 

Hi < — > H2 (11) 

The element 

f3=tlt2 (12) 

generates subgroup H3 of H. It is important that H3 is a normal subgroup of S, that is 
9 = ^^3 for any element g of S. Thus, we have the diagram of subgroups inside the 
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symmetry group S 

Hi 

\ 

Hs ^ n ^ S (13) 

/ 

in which the arrows signify the imbedding of subgroups. 

The group of symmetries, iS, plays the key role in finding minima of the potential V . The 
following general arguments, based on the theory of groups, arc quite useful in this respect. 
Consider a point // of space X of the angles 0i, 02, Suppose that is a minimum of U . 
Then points 

called the orbit of the point ^ under the action of the group S. arc also minima of U . 
The number of points [s! of the orbit can vary. In fact, let us consider all transformations 
g oi S that leave /i invariant, that is /i' — g ■ /i — /i. It is alleged to be known that 
the transformations form a subgroup of S, called stationary subgroup Hfj,. The stationary 
subgroups, TC^ and , for points fi and u of an orbit, are conjugate, that is = g~^ Hy g 
for an element g of S. The number of different points ^' equals to the ratio of the orders of 
S and 7Y^, that is to 2 or 4, depending on the choice of point /i. To be specific, consider a 
point /X having a stationary subgroup that coincides with the subgroup Ti. The latter is 
a normal subgroup of S of index 2, that is the factor set S/H consists of two elements. Thus, 
the orbit of /i under the action of S consists of only two points that correspond to the same 
value of U and have the same stationary subgroup Tl, because the latter is a normal subgroup 
of S. The situation is quite different if we take a point u having stationary subgroup Tii, 
which is different from 7^2- The subgroups do not coincide in S, even though they are 
conjugate. The orbit of /i under the action of S indicated above consists of four points that 
we may sort out as follows: two points having the stationary subgroup Hi and two points 
having 7^2- This is due to the fact that for one thing the subgroup Ti. is commutative and 
therefore its elements generate points of the orbit but with the same stationary subgroup, 
that is Til, and for another there is an element g that gives points of the orbit having the 
stationary subgroup 7i2- In contrast, a point having the stationary subgroup Tis has the 
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orbit consisting of four points which have the same stationary subgroup Hs, because the 
latter is a normal subgroup of S. 

IV. MINIMIZATION OF POTENTIAL FOR PAIR-INTERACTION. 

It is to be noted that numerical evaluation of the minima runs across a poor convergence 
of standard algorithms for minimization, because of rather flat surfaces of constant value 
for the function of three variables, U{4>i, 4>2, 0- To some extent, one may get round the 
difficulty by observing that for points that remain fixed with respect to a subgroup Q of 
S, the minimization problem is reduced to that for a smaller number of variables. This is 
due to the fact that for degrees of freedom perpendicular to the set of invariant points the 
necessary conditions for extremum are verified automatically and one needs only to study 
the conditions for longitudinal variables, that is to solve a smaller system of equations. 

To see the point let us consider a function f{x, y, z) of variables x, y, z even in x, so that 
f{x,y,z) = f{—x,y,z). The set of invariant points is y — z plane, and we may look for 
minima of the function f{x = 0, y, z), thus we need to solve only the two equations 

^f(x = 0, y, z) = 0, -^f(x = 0, y, = 

For larger groups of symmetries the number of variables necessary for calculations can be 
reduced even further, as happens for the minimization of C/(0i, 02, we are studying. It is 
easy to convince oneself that the sets of fixed points (0i, 02, 0; ^^^^ invariant under the 
action of a subgroup of <S, read as follows 

^1 ■■ (01 = 0, 02 = + TT, 
^2 : (01 = 0, 02 = - 0, 

^3 : ('/'i = ±f' '^2 = ±^, 
in which the J^i are invariant under the action of subgroups Hi, Ti.2, Hs, respectively. 

The analysis of symmetries of U given above enables us to sort out the minima according to 
the effective value of the phosphate charge Q, determined by the Debye-Hiickel screening. 
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FIG. 1: Charge Q against energy E for minima of U. Curves I, II, III correspond to minima of the 
types I, II, III, respectively. 




-7C 71 

FIG. 2: Cube of the symmetries indicating sets in space (f)i,(f)2,(, invariant with respect to subgroups 
of S. Main diagonal plane, B, corresponding to subgroup H2; two rectangles A perpendicular to 
B to Hi; solid lines 71, 74, and 72, 73 corresponding, to Tis and H, respectfully. 



It should be noted that minima of the pair-interaction U depend on the distance between 
molecules k, and the effective phosphate charge Q. The latter is the control parameter we 
employ in numerical simulation. It is also useful for the description of possible experimental 
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FIG. 3: Cube of the symmetries view from the top. Dotted Une describes invariant points 
corresponding to subgroup H^'-, dashed hne points !F\ and subgroup T^i; sohd circles 71 and 74 to 
subgroup 7^3 ; 72 and 73 to subgroup 7Y. 
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FIG. 4: Charge Q against ^ for minima of U. Curves II, III correspond to minima II, III. 

results. In this paper we are considering k to within 10.2 — 34 A. Effective charge Q of 
phosphate groups determines the neutralization; it varies to within — .6, in dimensionless 
units, Q = corresponding to the total neutralization. Charges Q that correspond to the 
charge inversion, have not been considered. The Debye length. A, has been varied to within 
7 — 35 A, depending on the ion strength of solution. 

The numerical data and the symmetry analysis given above suggest that there should be 
the following three types of the minima. 
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1. Type I characterized by the molecules having a cross-like conformation, " snowflakes", 
that is ^ being close to 7r/2. It exists for Q large enough. Its symmetry subgroup 
depends on the value of Q and may take values Hi, 0.2, Ti-s, H. Therefore, we may 
claim that there exist four sub-types of minima 1: It^^, It^j, 17^3, 1^^, each of them 
consisting of two subtypes which are given by specific conformations of the angle 
variables. 

2. Type II for which taking values to within 0.1°. Their symmetry subgroups are Tii 
and H2, each type consisting of two sub-types. 

3. Type III for which ^ is to within 1°, that is larger than for II. The symmetry subgroup 
is 7^3, and there are four constituent types of the same symmetry. 



V. LIQUID CRYSTALLINE PHASES OF DNA. 

The minima of the pair-interaction could provide a useful clue to the problem of finding 
possible conformations of large ensembles of molecules, even though the pair-interaction we 
have employed is a rough approximation for studying closely packed suspensions, it is as 
much as to say that the molecules in dispersions be packed in such a way that the pair- 
interaction should take on the lowest possible values. Following these guidelines we feel that 
the study of its minima could allow certain conjectures concerning liquid crystalline phases 
of the DNA. In particular, we may suggest that there is a number of different phases, some 
of them being the cholesteric ones. In this respect our results are in agreement: 

1. with papers 6j and 8| in which it has been verified that the cholesteric phases are 
related to small or almost zero effective charges of the DNA, 



2. with the theoretical conclusions of papers [12] - jlSj concerning the electrostatic mech- 
anism of generating cholesteric phases of the DNA, 



3. with paper [8| as regards the existence of polymorphism or multistability. 



We suggest that the description of cholesteric phases of the DNA should take into account 
the symmetries described by group S for the electrostatic interaction for pairs of molecules 
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of the DNA. Our results are described by a sequence of the embedded subgroups given by 
diagram f|T3|) . 



The symmetry classification of the minima, and possible phases, enables us to see differ- 
ences between liquid crystalline phases which at first sight appear to be physically identical. 
Thus the snowflake phase, which seems unique, if superficially estimated with the value of 
angle ^, may have several different types. The same is true for "cholesteric" phase II. In 
contrast, cholesteric phase III is unique, in the sense that its symmetry subgroup is only 
7^3. The subgroups provide a concise means for describing conformations of pairs of DNA 
molecules that determine the liquid crystalline phases. 

The physical parameter we use in our numerical work, is the effective charge Q, due to 
the phosphate charge of the DNA and the Debye-Hiickel screening. It enables us, at least 
qualitatively, to take into account the experimental setting. In fact, the charge distribution, 
mimicked by Q, can be effected by changing the constituency of solvent. The numerical 
results illustrated in FIG{T]are very suggestive. 

We see that for Q larger than a certain threshold value Qc the energy of the minima 
corresponding to phase I, called the snowflake phase, is smaller than that corresponding to 
phases II and III. Therefore we may suggest that for solvents that provide sufficiently high 
value of the effective charge Q there are no cholesteric phases proper and the snowflake phase 
could be present. 

It is surprising that liquid crystals with the symmetry properties described in this paper 
have so far remained unobserved. At least in part it might be due to the fact that the 
snowflake phase could be misinterpreted in the X-ray experiments. It is worth noting that 



2l|, in that it 



the formation of the snowflake phase is similar to that of the Wigner crystal, 
is due to the long-ranged Coulomb interaction of particles at low density, the inter-particle 
repulsion leading to the conformational organization. 

The phase diagram illustrated in FIG H] describes the behavior of the effective charge at 
constant inter-helical distance. Therefore, special precautions are to be taken so as not to 
change the latter while modifying the effective charge, with a view to see the transition of 
the snowflake phase into the cholesteric one. 
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It is also worthwhile to look for different cholesteric phases at small values of Q, that is 
less than Qc- The guideline being given by the fact that the size of differs by orders of 
magnitude for the phases II and III, at the same time values of the effective charge Q being 
equal. According to FIGS] one may expect phase transitions in which the cholesteric angle 
^ changes both its order of magnitude and sign. 

The formation of liquid crystalline phases of the DNA involves the degeneracy of S- 
symmetry, in accord with the general prescription of condensed matter theory, larger sub- 
groups corresponding to more isotropic phases. For the phases under consideration in the 
present paper it could be effected owing to the increment of frustration of the molecules 
trying to satisfy both the packing and electrostatic constraints. The symmetry considera- 
tions playing the key part, one may expect there could be phase transitions accompanied by 
changes in symmetry described by subgroups of group S. 

One may attempt to employ the conformation of charges of the molecule of the DNA that 
relies on the pyrophosphate groups rather than on the phosphate ones as usual. By now 
there has been a considerable progress in studying molecules of the DNA with substituted 



pyrophosphate internucleotide groups, j22|]. The samples studied in paper [22| had the double 
charges, —2 electron charge, due to the pyrophosphate groups inserted periodically, usually 
one group for every ten base-pairs of the DNA. One can expect that the study of dispersions 
of the DNA based on molecules carrying pyrophosphate internucleotide groups, could make 
for understanding the part played by electrostatic forces in generating liquid crystalline 



phases of the DNA. The design of novel DNA probes, [23|, provides new opportunities in 
this respect. 

The authors are thankful to S.Ya.Ishenko, S.A.Kuznetsova, and Yu.M.Yevdokimov for 
discussions. The work has been partially financed by the MSU Interdisciplinary Project for 
Studying Macromolecules. 
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